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Abstract 

We consider solutions to the full (non-isentropic) two-dimensional Eu- 
ler equations that are constant in time and along rays emanating from the 
origin. We prove that for a polytropic equation of state, entropy admis- 
sible solutions in L°° with non-vanishing velocity, density, and internal 
energy must be BV . Moreover, we obtain some results concerning the 
structure of such solutions. 



1 Introduction 

The study of multi-dimensional compressible gas flow (and more generally, 
multi-dimensional hyperbolic conservation laws) is a very active field of research, 
and there are many open questions. As a modest step toward understanding 
multi-dimensional flow, several reductions that are suggested by experiments 
can be made. In flow patterns such as regular and Mach reflection, the flow 
can be simplified to be constant in time and along rays emanating from a dis- 
tinguished point. A regu lar reflec tion consists of four shock waves meeting at a 
point, (see [CTTOl IEL081 ICKLOOl IZhe06l IMTC100I lElTTOl lElT09al IEU09 bj ) . while 



a Mach reflection consists of three shocks and a contact discontinuity meeting at 
a point (see [BD921 IBD061 IHor861 IHT021 IVK991 ISke97] ) . Therefore, we consider 
solutions U(t,x,y) that are in fact of the form U(6), where 9 = Z(x,y). Our 
goal is to investigate the regularity and structure of such solutions. 

Historically, understanding the Riemann problem was instrumental in the 
development of the theory of one dimensional conservation laws. In that case, 
one seeks solutions to the Cauchy problem in which the initial data is constant 
along rays emanating from the origin, which in one space dimension is of the 
form 
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Riemann problems are the foundation of theoretical existence and uniqueness 
proofs for general Cauchy problems (the Glimm scheme |Gli65| being the main 
example) , as well as of numerical schemes that accurately capture shocks (such 
as the Gudonov and Roe methods). 

The two-dimensional Riemann problem is much more complicated — some 
work has been done in which the initial data is constant in each quadrant (see 
|LZY98[ IZZ501 IZheOll IEE98] ). Instead of focusing on that kind of initial value 
problem, we seek steady in time solutions with more general configurations. 

Solutions of this form are also motivated by questions of non-uniqueness. In 
|E1106| , Elling numerically constructed an unsteady solution in which the initial 
data is a steady solution that is constant along rays. Therefore, understanding 
this class of solutions may lead to more examples of non- uniqueness, perhaps 
even analytical proofs. 

Another interesting issue for multi-dimensional conservation laws is finding 
the best function space to utilize. Whereas BV is the natural choice for one 
spatial dimension (see |Gli651 IGL701 IBBOlj ) , it is well known that it is not suit- 
able for the multi-dimensional case (see |Rau86] V However, when restricted to 
steady and self-similar solutions, the two-dimensional Euler equations are simi- 
lar to one-dimensional self-similar conservation laws, and so perhaps BV is an 
appropriate choice in this setting. This is one of the major results of this work 

- that under reasonable assumptions steady and self-similar Euler flows that 
are only assumed to be L°° can be shown to be BV. 

This paper is an extension of the work in |ER12j . but there are several 
important differences. In that case, the solutions were assumed to be small 
perturbations of a constant supersonic state, and so the steady problem was 
still strictly hyperbolic. In this case, we permit large variations in the solution. 
When the flow is sonic, there is a parabolic degeneracy in the steady problem; 
for subsonic flow the system becomes hyperbolic-elliptic. Therefore, a different 
approach from the one in |ER12j must be utilized. 

Moreover, the analysis in [ER12] required that the eigenvalues be distinct 

— in full Euler flow there is a double eigenvalue corresponding to shear waves 
and entropy jumps. Fortunately, many of the difficulties in dealing with this 
degeneracy can be avoided since when dealing with the Euler system we have 
explicit expressions for the shock curves and averaged matrices, and do not need 
to rely on the implicit function theorem. 

Sections 1-6 detail the various assumptions, derive the appropriate weak form 
of the equations, recall important properties of shock transitions, and demon- 
strate how to work with the solutions in the L°° setting. Section 7 shows that 
there can be a finite number of contact discontinuities, proves that shocks have 
neighborhoods on either side on which the flow is constant, and demonstrates 
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what can go wrong if the velocity vanishes. All that is required for these results 
is a standard equation of state. Sections 8 and 9 go into more detail about the 
configuration of shocks and prove BV regularity for the case of a polytropic 
gas. Finally, sections 10 and 11 describe results regarding the structure of pos- 
sible solutions in the polytropic case — how many waves can there be, what 
configurations can occur, etc. 

2 Preliminaries 

Let V <G IR" 1 be an open set. Consider smooth functions r), ip x , ip y ■ V — > IR. For 
A C IR 3 we say U = (t/\ U m ) € Ll oc (\R 3 ,V) is a weak solution of 

ri(U) t + V(U) X + ip y {U) y > in A, (2) 

if the inequality is satisfied in the weak sense: every x € A has an open neigh- 
borhood N so that for nonncgative smooth <f> with supp $ <e N, 

- [ $ t r){U) + ^ x i' x (U) + %^{U)d{x, y, t) > 0. (3) 

JlR 3 

We call U a strong solution (or classical solution) if, in addition, it is almost ev- 
erywhere equal to a Lipschitz continuous function. Weak solutions for a system 
of conservation laws 

Ut + .f x (U) x + f v {U) y = in A, (4) 

with f x , p : V — > IR m smooth, are defined by interpreting it as 2m inequalities 
with = replaced by < or >, and with r)(U) := U a ,^ x {U) = f xa {U),ip v {U) := 
f ya (U) for a = 1, ...m. We call (77, ip x , ip v ) an entropy-entropy flux pair for (01 
if 

1%=Wjf5, ru=Vuft on P. (5) 

A weak solution U of the Euler equations ((4]) is called an entropy solution if it 
satisfies for the physical concave entropy 77, which will sometimes be denoted 
by pS. The Euler equations are given by 





( p \ 




( P u \ 




( pv \ 


u = 


pu 


J x (u) = 


pu 2 + p 


J y (U) = 


puv 




pv 


puv 


pv 2 + p 




K p E ) 




{ u(E + p) J 




{ v(E+p) J 



where p is the density, u and v are the horizontal and vertical velocities, re- 

2,2 

spectively, p is the pressure, and E := p{ - ^ — h e) is the total energy per unit 
volume, with e the internal energy per unit mass. 

We now derive the weak form for steady, self-similar solutions U which only 
depend on 9, the polar angle Z(x,y). Integrating by parts in t eliminates the 
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first term in the integrand of (J5J), after which using compact-in- i support and 
integrating with respect to t yields the equivalent statement 

^ x iP x {U) + %ijj v {U)d{x,y)>Q (7) 

IR 2 

for all nonnegative smooth compactly supported t-independent functions $ : 
IR 2 — > IR. Now, change variables to polar coordinates to obtain 

OO />27T 

Jo 



0<- I I ($ x {rcos9 1 rsm9)ij x (U(9)) + $ y {rcos9,rsm9)il; v (U(9)) S jrd9dr. 

(8) 



Define a smooth 2tt periodic function </> as 

poo pQO pOO 

0(6*):= / <f>(r cos 9, r sin 6)dr = - cos 9 I $ x rdr — sin 6 l $ y rdr. (9) 
Jo Jo Jo 



Notice that 



/>oo />oo 

<t>'{6) = -sin0 I <S> x rdr + cos 9 / <5> y rdr. (10) 
Jo Jo 

We then have that 

POO 

<$> x rdr = cos 00(0) + sin 00' (0), (11) 



- / § y rdr = sin 00(0) - cos00'(0). (12) 
Jo 

Therefore ([8j becomes 

;00 + sin0 ( /)')i/. :l; ([/(0)) + fsin00-cos00')^([/(0))(20. (13) 



< / ( cos* 



This is the weak form of 

- sin 9ip x {U {9)) + cos 9ij.> v (U (9))) + cos 6ip x (U (9)) + sin 9ip y {U(0)) > 0, 

(14) 

or equivalcntly, 

sin 6tp x {U{6)) - cos0V y (^(0))) e - cos 6ip x (U (0)) - sin0^(£/(0)) < 0. (15) 

The left side is a nonpositive distribution, and is therefore a nonpositive Radon 
measure. By the Riesz representation theorem (see for example jBC09j . Chapter 
7), it is the distributional derivative of a non- increasing right continuous function 
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of bounded variation. Therefore, there is a version (that is, an element of the 
same L°° equivalence class) of U such that 



(sin 9ip x (U(6)) - cos9^ y {U(9)) 



< 



cos r]ip x (U (ij)) + smnip y (U(rj))dr). 

(16) 

Interpreting the conservation law as 2m inequalities, we obtain 

sm9f x {U{9)) - cosO f v (U {6) j) = cos9f x {U(9)) + sm9f y {U(9)). (17) 



The right side is L°°, and so the quantity being differentiated in the distribu- 
tional sense on the left must have a version that is Lipschitz. Therefore, the 
fundamental theorem of calculus holds and wc have 



sm0/ x (t/(0)) - cos9f y (U{9)) 



cos vf x (U(v)) +emrif v {U(ri))dt). 

(18) 



The following lemma from |ER12] shows there is a version of U such that these 
conditions will hold everywhere, not just almost everywhere. From now on, 
consider this version of U. 



Lemma 1. [ER101 Let f2 C IR™ measurable nonempty, K C IR m compact, 
U € L°°(fi) so that U{x) € K for a.e. x € 0, g : K -> IR* and g : Q -> \R k 
continuous. If 



g{U{x)) < g(x) for a.e. x e f2, 



(19) 



(meaning gi(U(x)) < gijx) for alii, where g = (51, ...,gk), g = (h, -,9k)), then 
we can find a version U of U , with values in K everywhere, so that 

g(U) < g for all x e Q. (20) 

3 Jump Conditions 

We have, for all 9\ and 6*2, 

-co S 9f y (U(9))+ fi m9f x (U(9))) 02 = 0{\6 2 - (21) 
On sequences {9^} — > 9 with U{9 ± ) — s- U±, we have 

sin 9 \f x (U)] - cos 9 \f y (U) } = (22) 
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(where [<?(?/)] = g{U+) — g(U-) for any function g of U). For the case of non- 
isentropic Euler, this becomes 



sin 8 [pu] — cos 0[pv] = 
sin 9[pu 2 + p] — cos 8[puv] = 
sin 8[puv] — cos 8[pv 2 + p] = 
sin 8 [u(E + p)] - cos 8[v{E+ p)} = 



(23) 
(24) 
(25) 
(26) 



where E 



-j^ — h ej is the total energy per unit volume, with e the internal 
energy per unit mass. 

To separate the cases of shocks and contacts, we introduce the normal (angular) 
and tangential (radial) velocities at 8: 



N := u sin 8 — v cos 8 
u = N sin 8 + L cos 8 



L := u cos 8 + v sin 8 
v = —N cos 8 + L sin ( 



We immediately observe that (|2"3"|) is equivalent to 

[pN] = 0. 

(|2"4")l is equivalent to 

sin6*[p] — [cos 8 puv ~ sin^pu 2 ] = — [puN]. 
Similarly, j25]) yields 

cos6*[p] = [sin 8 puv — cos6*pu 2 ] = [pvN]. 

Therefore, 

= sin 6[pvN] + cos 8[puN] = [pNL]. 



(27) 
(28) 



(29) 
(30) 
(31) 
(32) 



This means that p+N + L + - p_N-L_ = 0, or p + N + (L + - L_) = (from ([29]) ). 
Therefore, if p+N+ = p-N- ^ 0, [L] = 0. Hence, if there is mass flux through 
a shock, the tangential velocity is continuous. However, if N + = N~ = 0, then 
the tangential velocity may be discontinuous (called a contact discontinuity). 
Also, we have 



[p] = - sin 8[puN] + cos 6[pvN] = - [pN 2 ] 
Finally, (I26|) is equivalent to 

"1 



[pN 2 + p] = 0. (33) 



0= [N{E + p)\ 



pN\u\ 2 + pNe + Np 



For the case of a shock, divide by p+N + = p_iV_ to obtain 

T 



-\u\ +e+pr 



= 0, 



(34) 



(35) 



G 



where r := p 1 is the specific volume (volume per unit mass). Denote p_iV_ 
p+N + =: M. Then j33]) becomes 



]p] = -M[N]. (36) 

Therefore, 

-M^-Mt^^-MK (38, 

P-P+ 

Multiplying (06} by (r_ + r+) yields 

[p](T_ + r+) = -M[N](t- + r+) (39) 
= M(N- -N+) P ~ +P+ (40) 



M(p-N_ - p + N+ - p-N + + p+N-) 



(41) 



P-P+ 

M{t-N- - t+N + ) = Nl - N'l = -{\u\ 2 }. (42) 



Substituting this into ([33)) yields 

[e + pr] = -\p](T-+T+), (43) 

[e] = ~(p_+p+)[r]. (44) 

Therefore, for a shock wave with fixed (r+,p+), the states (r_,p_) that can 
be connected by a shock are defined by H(t~,P-) = 0, where H(r,p) is the 
Hugoniot function defined below: 

H(r,p) = e(r,p) - e(r + ,p + ) + -(r - r + )(p + p + ). (45) 



4 Weak form in terms of tangential and normal 
velocities 

Analogous steps to those used to simplify the jump conditions to be in terms 
of radial and angular velocities can be done on the weak form of the equations. 
Multiplying a distribution (in this case an L°° function) by a smooth function 
(sin 9 or cos 9) results in another distribution, and the product rule applies for 
distributional derivatives of distributions multiplied with smooth functions. 
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(p~7|) with the Euler fluxes becomes the following system: 

sin 9pu — cos 9pv J = cos 9pu + sin 9pv, (46) 



(^sm 9(pu 2 + p) — cos 9 (puvfj = cos 9 (pu 2 + p) + sm 8(puv), (47) 

^sind(puv) — cos9(pv 2 + p)^j — cos9(puv) + sm8(pv 2 +p), (48) 

(sm6(u(E+p)) - cos 8 (v(E +p))\ = cos9(u(E + p)) + sm8(v{E +p)). 

(49) 

Substituting the definitions of N and L we obtain: 

(pN) g = pL, (50) 

(upN + psai6) g = upL +pcos9, (51) 

(vpN — pcos9)) e = vpL +ps'm9, (52) 

(N(E+p)) 8 = L(E+p). (53) 

Note that 



sin 9 [upN + psin6*) e = (us\n9{pN) + psin 2 9) e — cos 9{upN + psin6>), 
cos 9{ypN — pcos #)) _ = ( — v cos 9(pN) + pcos 2 9) s — sin6 (vpN — pcos 9)) . 



(54) 
(55) 

Thus, sinl? • (J5TJ - cosl? • (J52]) yields 

(pN 2 +p) g - P NL = pNL, (56) 

or 

{pN 2 +p) g = 2 P NL. (57) 

Similarly, 

cos 8 (upN + psuv8) g = (u cos 8{pN) + p sin 9 cos 9) + sin 8 (upN + p sin 8) , 

(58) 

sm8(vpN — pcos 6*)) e = (w sin8(pN) — p sin 8 cos 8) „ — cos 9{ypN — pcos 9)) , 

(59) 

and so cos 9 ■ d5T]) + sin 9 ■ ([52]) yields 

(pLN) g + pN 2 + p = pL 2 + p, (60) 



or 



(pLN)g = pL 2 - pN 2 . (61) 
Therefore, the Euler equations are equivalent to 

(pN) e = pL, (62) 

(pN 2 +p) e = 2 P NL, (63) 

(pLN)g = pL 2 - pN 2 , (64) 

(N(E+p)) e =L(E+p), (65) 

satisfied in the distributional sense. 

5 Shock sides 

For now, consider self-similar solutions U = (p, pu, pv, E) : S 1 — > V £ IR 4 that 
are L°°, where the phase space V consists of states U satisfying: 

< C _x < p < C < oo (66) 

< C- 1 < e < C < oo (67) 

x= < C < oo (68) 



where C is a positive constant. Recall E := p ^ - ^ ' — h ej , where e is the specific 
internal energy. Denote by p the pressure, T the temperature, r the specific 
volume, p = - the density, and S the specific entropy. From thermodynamics, 
knowing any two of p, T, r, e, S 1 determines the other three. 
We assume that the equation of state for pressure is given by p = g(r, S) and 
that it is smooth for the phase space under consideration. We also assume 

g T := -p 2 c 2 < 0, (69) 
9rr > 0, (70) 
9S > 0, (71) 
where c > is the sound speed. Furthermore, from thermodynamics we have 

de 



dr 
de 
dS 



= -P (72) 

S— const 

= T, (73) 



or equivalcntly 



(t) - ~ 2 (74) 



Of 

as 



= T. (75) 

p— const 
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To emphasize that this definition of sound speed, 



®P J S=const 

agrees with another definition seen in the literature, 



p f dp 



dp 



p— const 



consider 



(76) 



(77) 



p(p, S) =p(p,e). 



(78) 



Differentiating with respect to p while holding S constant, we obtain 



dp 

dp 



S— const 



dp 



dp 
dp 



S— const 



dp 
de 



p— const 



de 
dp 



S— const 



dp 
de 



(79) 
(80) 



p— const 



We now argue that even in the L°° setting, in which left and right limits may 
not exist, there still exists a well defined notion of shocks and contacts, as well 
as front and back sides of shocks. For the remainder of this paper, a subscript 
max or min on a quantity refers to the maximum or minimum permissible value 
of that quantity for U € V. 



Lemma 2. If U is discontinuous at 9 , then it can be well defined as having 
either a forward facing shock, a backward facing shock, or a contact discontinuity 
at 0q. 

Proof. Suppose U (9) is discontinuous at 9q. Then we can pick sequences {9 n } , {6' n } 
8 with U{9 n ) — > U , U{9' n ) — > U' 0l with U ^ U' Q . Then, from above, necessarily 
PqNq = PqNq. Based on our assumptions about p, it follows that N and iVg 
are either both positive, both negative, or both zero. 



Case 1: Nq = Nq = 0. Let {#"} be any other sequence converging to 9o, and take 
any subsequence 0'^ k y Since U is L°° , there exists a subsequence 
such that U(6'^ k ,j\\) — > Uq . Applying the jump conditions to Uq and Uq shows 
that ./Vp' = 0, because -/V = 0. Therefore, since any subsequence has a subse- 
quence converging to zero, N{9' 7 [) — > 0, and as {#''} was arbitrary, we see that 
in fact N(9) is continuous at 9q and N(9q) = 0. In this case we say that U has 
a contact discontinuity at 9q- 
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Case 2: No,Nq > 0. We claim that there exists some neighborhood contain- 
ing #o on which N(0) is positive. Suppose not. Then there exists — > 8 

such that N{0' 1 ' 1 ) < for all n. Again, there exists a subsequence such 
that U(6'^, k s) — > Uq, and by assumption Nq < 0. However, applying the jump 
conditions to No and Nq gives a contradiction. Therefore, N(0) > on some 
neighborhood of O , and we call the shock forward facing. (The interpretation 
is that gas particles enter the front side of the shock, and leave the back side. 
Since the flow of mass through the shock is aligned with our choice of normal 
vector for the shock, we call this case forward facing.) 



Case 3, for which Nq,N± < 0, is similar. Wc call this case backward facing. □ 



6 Entropy 

The entropy admissibility criterion, with 9~ < 0+, takes the form 

(- cos6^ (U(9)) + sm6r(U(6)))\ e s i+G(\6+- 6- \)<0. (81) 

(Note this is reversed from the usual conservation laws literature since we are 
using a concave (physical), not convex (mathematical) entropy.) In the limit 
0% -> 8, if U(0n) -> U±, this becomes 

sm0[^ x (U)} - cos9[ip y {U)} < (82) 

For the case of the Euler equations, this reads 

> sin 0[puS] - cos 6[pvS] = [pNS], (83) 

which becomes 



M[S}<0. (84) 

(where Ai = p_N_ = p + N + ). In the case of a forward facing shock, S + < 5L, 
and in the case of a backward facing shock, 5+ > S-. In either case, the entropy 
cannot decrease when the gas passes through the shock from the front to the 
back. 

We now assume the case of a forward facing shock, and fix a pair of left and 
right sequences 6„ as usual, so that 9~ < 6*+ for all n, and U(0^) — > U±. Wc 
know that U- must satisfy H(t-,P-) = for the Hugoniot function for state 
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Lemma 3. ( \CF48j ) Assume that for fixed (r + ,p + ), the set of states (r,p) for 
which H(r,p) = is a smooth curve in the (r,p) plane that can be described 
by p = G(t). Then, for the case of an entropy admissible forward facing shock, 
which requires S- > S+, we have 

1. r_ < r + (that is, shocks are compressive), 

2. N+ > c + > 0,0 < iV_ < c_. 

We have the analogous statement for an entropy admissible backward facing 
shock, which instead requires SL < S+: 

1. T_ > T+, 

2. -N- > c_ > 0,0 < -N + < c+. 
Moreover, for any shock, 

\N±\ > CmmPmm . (85) 

Praax 

Proof. This proof can be found in |CF48j . First we derive some properties from 
the equation of state p = g(r, S). Recall 

g T = -c 2 p 2 < 0, (86) 

9rr > 0, (87) 

g s > 0. (88) 
Differentiating S = S(r,g(T, S)) with respect to S we obtain 

1 = S p g s => S p >0, (89) 
from (|55|) . Differentiating S = S(t, g(r, S)) with respect to r yields 

= S T + S p g T S T > 0, g T = (90) 

Dp 

from (JHH1) an d ([56]) . Differentiating again with respect to r yields 

= S TT + 2S Tp g T + S p g TT + S pp g 2 T . (91) 
From ([87]). it follows that 

S TT + 2S Tp g T + S pp g T < 0. (92) 

S 

Substituting in g T = — — and multiplying by S p yields 
o p 

S TT Sp *2S Tp S r S p ~t~ S pp S T <^ 0. (93) 
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We now claim that G"(r + ) > 0. We start differentiating the Hugoniot function 
with respect to t along the curve p = G(r). 

e(r, S(t)) - e(T + ,S(r + )) + \{t - r + )(G(r) + G(r+)) = (94) 

e s S'( T ) +e T + 1 -{G(t) + G(r + )) + i( T - r+)G'(r) = 0. (95) 
Substituting in es = T, e T = — p yields 

rS'(r) + i(G(r+) - G(r)) + \{r r + )G'(r) = 0. (96) 
Therefore, S"(r + ) = 0. Furthermore, 

T'{t)S'{t) + TS"(r) - \g'{t) + \g'{t) + \{t r + )G"(r) = 0, (97) 

implies S"(t + ) = 0. The relation G(r) = ,g(r,5(r)) yields 

G'(T)=g T +gsS'(T), (98) 
G"(r) = 5rT + 2g TS S'(T) + gssS'ir) 2 + g s S"(r), (99) 

and so G"(t + ) = g TT (r + , S(t + )) > from our assumptions about the equation 
of state. 

Now, consider any other state (r_,p_),p_ = G(r_), such that 7J(t_,p_) = 0. 
Parametrize the straight line segment between these two states as 

T = T+ + (T_ - 7+)s, P = _p+ + (p_ - _P+)S. (100) 

(H(t) and 5(r) refer to states on the Hugoniot curve, while -ff(s) and 5(s) refer 
to states on this straight line segment.) We have that 

H'(s) = TS'(s) + \({p + -p_)(r_ - t+) - (t_ - r + )(p_ = T5'(a). 

(101) 

Hence if H'(s) = 0, S'(s) = 0. i?(0) = H(l) = 0, and so H'(s ) = for at least 
one So between and 1. Therefore, S'(sq) = for at least one So between and 
1. However, consider S"(s). 



V- -p+ 



S'(s) = S T r'(s) + S pP '(s) ^ = - 7 7 - (102) 

<->D 



s=s t_ — r 4 

2 



S"'(«) = 5 TT r'(s) 2 + 2S tp t'(s)p' (s) + S pp p' (sf (103) 

= (P- - P+) 2 (Srr { ^~ T+) l + 2S Tp T -~ T+ + S pp ) , (104) 

V (p--p+) p--p+ J 



S"(s ) = (P J +)2 (S TT S 2 p - 2S Tp S T S p + S PP S 2 T ) < 0. (105) 
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Therefore, 5(s) has a single maximum at so and no other critical points. There- 
fore, 5'(0) > 0, 5'(1) < 0. Therefore, H'(l) < and this straight line cannot 
be tangent to the Hugoniot curve at s = 1, or cquivalently r = t_. 
Now, consider the Hugoniot curve. We have that 

F'(r_) - TS'(r-) + ^((P+ -P-) + (r- - 7+)Gf'(r_)) = 0. (106) 



Suppose S'(t-) = 0, then 



G'(t-) = P ' P+ - (107) 



That is, the Hugoniot curve is tangent to the straight line between (r+,p+) and 
(r_ ,p-). This contradicts what we showed above, and so 

S'(t) ^ for aU r ^ t+. (108) 

Differentiating (p?T)) we obtain 

T"{t)S'{t) + 2T'(t)S"(t) + TS"\t) + i(( T - r + )G'"(r) + G"(r)) = 0, 

(109) 

and at r+ this yields 

S'"(t+) < 0. (110) 

Thus for (r_ — t + ) sufficiently small, 

SL > 5+ r_ < r + (111) 

5_ < 5+ r- > r + (112) 

which allows us to conclude 5'(t + ) = 0, S'(t) < for all t/t + . In conclusion, 
for a forward facing shock so that r + is a limit from the front side, 5_ > 5 + if 
and only if r_ < r + , and we can further conclude that SL > 5+ if there is an 
actual discontinuity. 

Finally, since S s = S p (p- -p + ) + 5 T (T_ -t+), the conditions 5(0) > 0, 5(1) < 

5 

when combined with g T = — — - = —{?<? and S p > yield 

Op 

(p_-p + )-(p + c + ) 2 (r_-T + ) >0, (113) 
(p_-p + )-(p_ C _) 2 (r_-r + ) <0, (114) 

leading to (taking t_ < t + ) 

( P+C+ ) 2 < _^+^ < (p _ c _)2. (ii5) 
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J35J yields 



(p+c+f < m 2 < (/J_c_) 2 , 



(116) 



from which we obtain Lax-type conditions N + > c + > 0, < N- < c_. 
The case of a backward-facing shock is similar, and the remaining lower bound 
on |iV±| follows from the fact that one must be greater in magnitude than c m ; n , 
and the relation 

N ± =N T ^. (117) 

□ 



7 Properties of discontinuities, uniform distances 
between different types of discontinuities 

We now use the idea of averaged matrices, that is smooth matrix valued func- 
tions of i7 ± and 8 that satisfy: 

(a) smd(r(U+) - f x (U-)) - cos6(fy(U+) - P{U~)) 
= A(U-,U+;8)(U + -U~) 

(b) A(U,U;9)=A(U;9) := mn9f§(U) - cos 9f^(U), 

(c) A(U~ , U + ; 9) is diagonalizable with real eigenvalues for all ?7 ± € V . 9 G 
[0,2tt[. 

The existence of matrices with the important conservation and diagonaliz- 
ability properties is guaranteed (see |H83j ) if our system possesses an entropy- 
entropy flux (77, ij) x , ip y ) with rj positive definite, but one that is easier for com- 
putations in the polytropic setting uses a different line integral in phase space 
as opposed to one using the entropy gradient as described in |H83j . It has the 
property that the averaged matrix is simply the flux matrix evaluated at a suit- 
able averaged state. This is called the Roe linearization for the system and is 
very common in numerics. 
Rearranging (|21[) , we obtain 

sm9 1 ( K r(U(9 2 )) -f x {U(9 1 ))) -coB0i(/»(tf(0 a )) =O(|0 2 -0i|). 

(118) 

Using the averaged matrix, we have 

A{U{8 2 ),U{9 1 )-8 1 ){U(8 2 )-U{8 1 ))=0{\8 2 -8 1 \). (119) 
Now, consider the matrix 

A(U;8) = sm9ff}{U) - cos9f^(U). (120) 
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By direct computation, we have that 

det(sm6ffi(U) -cos0/$(tT)) =0 <=> \N\ = c or N = 0. (121) 
We then have the following theorem. 

Theorem 1. Suppose U is continuous on an interval ]Oi,02[ and that \N\ ^ ±c 
or on this interval. Then U is constant on this interval. 

Proof. Fix some 9 G]#i,02[- We claim that U must be Lipschitz at 9. Suppose 
not. Then we can choose a sequence {9 n } —> 9 (with \0 n — 9\ ^ 0) such that 



< 



U{9 n ) - U{9) 



Divide both sides of CDS]) by \U(9 n ) ~ U(9)\ to obtain 
iW,*);,)/*)-"") 

1 



\U{9 n )-U{9)\ 



\U{0 n )-U{0)\, 
0(\6 n -9\) =o(l) asrwoo. (122) 



By assumption, sin Of ft (8) — cos 9fj^(9) is regular, so for |0„ — 9\ sufficiently 
small, A(U(9), U(9 n ); 9) will be uniformly regular, because the eigenvalues of A 
are continuous functions of U (see for example jSlOj ). and U is continuous at 
8 by assumption. That is, 



3d > Vz G IR" 



i(C/(#), [/(#„); 0!))z| ><$|z| 



Taking n — > oo, the left hand side stays bounded away from zero, while the 
right hand side goes to zero, leading to a contradiction. 

Therefore, U must be Lipschitz on ]0i,02[- Assuming 9 is a point of differ- 
entiability of U, we obtain 

im0f$(U(ff)) - cos0 f*(U(0)))Ug = 0. 

However, as we assumed the matrix was regular on ]9\ 1 02[, it follows that Ug = 
on this interval. A Lipschitz function is the integral of its derivative, so U is 
constant on ]9i, 02 [• d 

Theorem 2. Suppose there is a shock at 9a. Then there exist a + (9o) > 0$ and 
a~ (9q) < 9q so that U is constant on ]a~(9o), 9q[, ]0q, cr + (9o)[- In particular U 
has well defined left and right limits at shocks. 

Proof. We consider the case of a forward facing shock; backward facing shocks 
can be treated similarly. If there is a shock, choose left and right sequences 
9f L — > 9 with U{9^) — > U± with 6~ < 6*+ for all n. We have from above that 

N + > c+ (123) 
c_ > > 0. (124) 
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Suppose there is no <j + (8 ) such that N > c for all 9 &]9 , a + (9o)[- Then pick 
a new sequence #+ + \ 9 such that (passing to subseqences if necessary) 

N(0++) < c(9++) (125) 
U(9++) -+ U++ (126) 
6+ < 9++ (127) 

for all n. If C/++ 7^ f/+, then the Lax condition requires TV++ > C++, which 
contradicts our construction. If £/++ = [/+, then -ZV++ = TV+ > c+ = C++, also 
a contradiction. Therefore there exists a + (9o) > 9q such that N > c for all 
9e}9 ,o-+(9 Q )[. 

The fact that there is some <j'~(9q) < Qq such that TV < c for all 9 £ 
]cr'~(#o): #o[ is analogous. 

As observed in the proof of Lemma [3J N > for some neighborhood of 
$o, otherwise the jump conditions cannot be satisfied. Taking er~ (9 ) to be 
the left endpoint of this neighborhood intersected with ]a'~(9 ), 9 [ shows that 
0< N < c for all e]o- - 0o,0o[- 

All that is left to show is that U cannot possess any other discontinuities in 
these left and right neighborhoods. Contacts are already ruled out - we know 
that N = at a contact. Therefore we must show there cannot be shocks. 
Suppose there was a shock at 61 £]#o, c + (#o)[- It too must be forward facing, 
since N is positive for all such 9. Then, for some r\ &]9 ,a + (9o)[r\}(j^(8 1 ),9i[ 
we would have 

N(r)) > c( V ) < N( v ), (128) 

a contradiction. Therefore, U is continuous on ]9o, (J+ (9o)[ and Theorem [TJ shows 
it is constant. Similarly we can conclude U is constant on ]a~(0o), 9q[. □ 

Similarly, contact discontinuities are isolated and must also possess constant 
neighborhoods. 

Theorem 3. Suppose U has a contact discontinuity at 9q. Then, either (a) 

there exist o~ + (9q) > 9q and o~~ (60) < 9q such that U is constant on ]o~~(9q), 9q[, ]9q, o~ + (9q)[, 

or (b) 9q is contained in a closed interval on which \u\ = 0. 

Proof. From previous analysis we know that in fact TV is continuous at 9q and 
that N(9q) = 0. Therefore, choose a ir > 6 > such that 

\9-9 Q \<5 =^ \N(9)\ < Cminftnin . (129) 
, 2p 

max 

By the last part of Lemma |3l we have the lower bound for normal velocity at a 
shock: 

|/V ± 1 > Cmi "^ mi " , (130) 

Pm&x 
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Therefore, for | ^ — 6*o | < 8 there can be no shocks. If there is another contact 
discontinuity, N must still be continuous, and so N is continuous for 1 6* — I < 
Therefore, the set 



C := 



\e-6 \ < 5,N(9)=0 



} 



(131) 



is a closed set in ]9o — S, 9q + 5[. Its complement in ]9q — 5, Oo + S[ is therefore 
open and is thus a countable union of open intervals. If one of these open 
intervals ]ai, bi[ has neither endpoint not equal to 9q ± S, then U is constant on 
(oi, bi), and so U = Ui on ]aj, bi[, and N(a,i) = N(bi) = 0. However, since N is 
continuous we can take limits and find 



Then, since ^ bi, the vectors (sin cij , — cos a.; ) and (sin bi, — cos 6i) span IR 2 , 
and so = 0. However this contradicts that TV 7^ on Ja^&J. There- 

fore, C is either a closed interval containing 9 or simply {6* } . Consider the 
conservation of mass equation, 



in the distributional sense. If TV = on a closed interval, we can take a strong 
derivative at any 9 in its interior to obtain that L = on the interior. Therefore 
the supposed contact at 9q could not have had a jump in tangential velocity, 
only density and entropy (continuity of pressure is required by the jump con- 
ditions). (This opens the possibility to very irregular solutions. If the velocity 
field is zero on some interval, then a completely arbitrary density distribution 
can be prescribed, as long as the entropy /internal energy /temperature is also 
prescribed to result in constant pressure). If we disregard this possibility, then 
it follows that this closed interval must be a single point, and so there can only 
be one contact for \9 — 6 \ < 5. These very irregular solutions are not all that 
surprising, considering that when the velocity is identically zero the Eulcr equa- 
tions become pe = in the sense of distributions, so p = constant. Note that 
for isentropic flow, this situation cannot occur, since constant pressure can only 
be attained if density is constant. □ 

Note we have shown that the set on which U has a discontinuity is countable 
and discrete. Therefore, right and left limits are well defined, and so from here 
on out we modify U to be right continuous at every point. 

Lemma 4. Assume there are no stagnation points (that is, \u\ ^ everywhere) . 
Then the set on which N = is a finite set of points. Moreover, if N(9q) = 0. 
and there is a shock at 9' , 



N(a,i) = Ui sinai 
N(bi) = Ui sinbi 



Vi cos a; = 
Vi cos bi = 



(132) 
(133) 



(pN) = P L, 



(134) 



|0'-0o|><5>O, 



(135) 



with S independent of U . 
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Proof. Take any 9 with N(0 ) = 0. Define 



J ) ■= sup ^ 77 

»7>0o 



< |TV(0)| < Cl " inPmi " \/6€}6 ,vl}- 



This suprcmum is defined because it is taken over a nonempty set by Theo- 
rem [3] Moreover, since there cannot be any shocks on ]6q, a + (6o)], it follows 
that \N(a+(d ))\ = or c 2^'° ■ We claim that if <r + {9 ) - 6 < n, then 
\N(a + (9o))\ = c *2p Pmi ° ■ If not, then since there could be no shocks or contacts 
between and <r + (8 ), by Theorem [TJ U would be constant. However, this 
would lead to a contradiction of N(8 ) = N(a + (9 )) = and they were sepa- 
rated by less than n. Therefore, \N(a + (9o))\ = c ^ pmin , and U must be constant 
on }9q,o- + (9q)} (and continuous at a + (9o)). When U is constant, — L, and 
so 

^P^ = \N(<j+(9o))-N(9o)\ (137) 

= lim / L{0)d0 <\u\ max (o-+(9 )- 9 ). (138) 

Therefore, (a + (9o)— #o) is bounded below independent of U, and another contact 
could only happen for 9 > a + (9o). A similar argument works for 9 < 9o, and so 
the total number of contacts is finite. 

This same calculation shows that the distance between a contact and a shock 
is lower bounded independent of U. □ 

Lemma 5. If there is a forward facing shock at 9, and a backward facing shock 
at 9' , then 

|0-0'|>tf>O (139) 

for some 8 independent ofU. 

Proof. For forward facing shocks, the normal velocity is positive on cither side. 
Similarly, it is negative on either side of a backward facing shock. Therefore, 
between a forward facing and backward facing shock, TV must transition through 
zero, not jump between positive and negative values. In the previous lemma, 
we showed that the distance between a point at which TV = and any kind 
of shock is uniformly bounded away from zero (independent of U), and so the 
claim follows. □ 

8 Shock strengths and neighborhood sizes 

We now turn to the case of a polytropic gas, for which 

e = ra (140) 
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We assume that 7 > 1, recalling for air 7 = 1.4. We now find expressions for 
various quantities at a shock, following }T09j . Substituting the expression for e 
into we obtain 



\pN{N 2 + L 2 ) + ^-Np 
2 7 — 1 



However, [pNL 2 ] = for a shock, and so we have 



0. 



1 



-pN J 



1 



■Np 



= 0. 



.2 r 7-I 
This can be rewritten as (recalling p~N- = p+N^ 

1 



= h p+ N + )Nl + ^L-N + p, 
2 7 — 1 



However, since 



we obtain 



7V_ = N+ + 



:(p+N + )N± 



P+ -P- 



7-1 



N_p_ 



p+N^ 



-- ^(p + N + )Nl + ^-N +P+ 
7-1 



(p+N+) N + + 



P+ - P- 
P+ N + 



7-1 



P-IN+ + 



P+ -P- 
P+N+ 



(141) 
(142) 
(143) 
(144) 

(145) 
(146) 



Then, 



7 



N+(p+ - p_) - iV + (p+ - - 



Hp+-p-) 2 



7-1 

P+ - P- 
p + N+ V7-1 



2 P+ iV + 



1 w2 !/ x 7 



7 P+-P- 

tP-- 



p+ -p- 

~ (7 - l) P+ N+ 

P+ -P- 
" (7 - l)p+N+ 

Therefore. 



P+N'i 



7-1 



7-1 
(P+ - 1P- 



2 7 — 1 7+1 
p+N, — p + — p. 



Nl = ^-[p+{ 1 -l)+p-{ 1 + l) 



2P4 

P± f 7-1 g- (7 + 1) 
P+ V 2 p+ 2 



7-1 p+iV+ 

(147) 

(148) 
(149) 
(150) 



(151) 
(152) 
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Introduce 

V- -p+ 



z = 



> (153) 



P+ 

as the shock strength of a forward facing shock, so that 

P-=P+(l + z). (154) 

The remainder of these calculations are for forward facing shocks, and the same 
statements hold for backward facing shocks but with ± switched, and all normal 
velocities are negative. 

T 2 _ P+ ( 7 - 1 . 7+1 

= — { l + z —) (156) 

+ (157) 
Therefore we have the relations (assuming a forward facing shock) 



^i = — [^— + ^—(l + z)\ (155) 



1 + z^ (158) 



7V + - C+ = c + ^l + z^i-lj. (159) 

We need similar relations for the state behind the shock. 

N- = (p+Nl -p-+ p+) (160) 

P+N+ 

1 



p+N + 

l + z ±-^-^\ (162) 



2 

ci 



iV. 
iV. 
7V + 



4 (i + 


7 + 
Z 2 7 


7+1 


zp+ 


27 




7 + 1 




27 




7 -l\ 




27 J 





-p + (l + 2 )+p+ (161) 



+ ' l + Z -LJ_1 _ ~ ) (163) 
+ 

' l + s-i—l). (164) 
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Next, we obtain expressions for p_, then c_, iV_ — c_, and — — . 



P- = 



p+n+ p+m 



i 



N- 



l + z 



7-1 
27 



P+ 



1 + ^ 



IP- _ 7P+(1 + 2) 1 + zl W 



P- 

1 + 

27 



, 7+1 

27 



1 + z 



7V_ 



1+ 



1 + 2 



7-1 
2 7 



(1+^) 



1 + z 



7-1 
27 



1 + ^ 



1 + 

2 7 



1 + *^ 



C+\ 

7+1 \ 

2-v ' 



(1+*) 



\ 



1+^ 



1 + Z^— ^ - Vl + Z 



C- 



- 1 = 

C- 



1 + ^ 



1 + *^ 



(165) 
(166) 
(167) 
(168) 

(169) 
(170) 

(171) 

(172) 
(173) 



Vi + z 

We now argue that z is a suitable measure of shock strength. 

Lemma 6. At a shock, [\U\] and z are equivalent measures of the strength of 
the shock, i.e. there exists C > such that 

J < [\U\] < Cz (174) 

Proof. First, recall that we assume pressure is bounded away from and 00, so 
we only need estimates valid for 

Pmax Pmin 



(|166[) shows that 



< Z < Zr> 



M\ 



Pn 



, 7+1 
' 2 7 



= p+- 



,7-1 • 



(175) 

(176) 
(177) 
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The function 



7 + 



is concave for z > 0, with derivative at z — equal to -. Therefore 



Pmir 



Clearly 



<IMI< 

Finally, since — m[7V] = [p], and 

PminCmin ^ |^| ^ Pmax | ^ | max 5 

£*min - i r a ri i ^ T^max 
— Z < \[N\\ < Z. 



^max ^ max 



(178) 

(179) 

(180) 

(181) 
(182) 



We have therefore shown that in terms of the primitive variables V := (p, u, v,p) 
that for a shock 



C 



<\[V]\<Cz 



(183) 



for some C > 0. U = (p, pu, pv, pE) has derivative 



/ 1 

u p 

v p 

v «!±2! pu pv _i_ y 



(184) 



Therefore it is C 1 , and the operator norm of its derivative is bounded uni- 
formly on the phase space under consideration. We can assume our phase space 
in V variables is convex (if not we can still obtain a Lipschitz estimate from 
Jq Uy {s(V + — VJ) + V-)ds(V + — VJ) since U(V) have derivatives with operator 
norm bounded away from and oo for < C~ l < p < C < oo) and therefore 
this is a Lipschitz map. Furthermore \[U]\ > \[p]\ > fj, and so putting it all 
together, there must exist C > such that 



C 



< \[U]\<Cz 



(185) 



for a shock. We will thus denote J(U;6) := \[U]\ as the size of the jump in U 
at 9. □ 

We now estimate the sizes of the neighborhoods on either side of a shock on 
which U must be constant. 
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Theorem 4. Suppose U has a forward facing shock at 9q ■ Then the o~ + (9q) > 9o 
and o~~(9q) < $o (from Theorem^ such that U is constant on 
]cr~ (6q),6o[,]9o,o- + (#o)[ satisfy the following: 

o- + (9 )>9 + 5 L J(U;9 ) (186) 
o--(6 Q ) <6o-5 l J(U;0 ), (187) 

where Sl is a positive constant independent of U . Furthermore, 

N+-c + >S L J(U;6 ), (188) 
N--c-<-5 L J{U\6 Q ). (189) 

The analogous statement holds for a backward facing shock. 
Proof. Suppose the shock is forward facing. Take 

o-+{9 ) := sup in N{9) - c{9) > V0 e]6> , //[} . (190) 

The supremum is taken over a non-empty set from Theorem [2j and is thus 
defined. Furthermore, it must be the case that N(a + (9o)) = c(9q) — from 
Theorems [2] and [3] there could not be a shock or contact at a + (9o), and so U is 
continuous at o- + (9q). Furthermore, it is constant on ]9q,o + {9q)[ by Theorem 
[TJ (Note however there could be another forward facing shock at 9' > a + (9 ) 
arbitrarily close to ct + (9q)). When U is constant, = i, and so we have 

rv 

N + - c+ = N(9 Q +) -N(a + (9 )) = lim / L(0)# (191) 

*7\t0 o J<7+(e ) 

<(a + {6 )- 9 )\u\ max . (192) 

However, recall 



N + c + = c + (J 1 + z2±i _ 1 j > Cmin + z l±± 1 j . (193) 
The function on the right is concave for z > 0, at z = 0, and therefore 



7+1 

iV+ - c+ > c min AV Z L z , (194) 



All together, we obtain 

(<7+(0 o )-0o) >6' L J(U;0o), (195) 
for some <5^ > 0. Similarly, define 

<r~(0 o ):= inf U < N{9) < c{9) V9 e]r),6 [\ ■ (196) 
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As before, this infimum is over a non-empty set by Theorem[5] Similar to before, 
it must be the case that either N(a~(9 )) = c(<r~(9 )), or N(a~ (6 )+) = 0, 
since there can be no shocks on [a~ (6o),0q[, and no contacts on ]a~(9o), 9q[. 
However, if N(a~(9 )+) = 0, the only possible discontinuity would be a contact, 
at which N is continuous, and so N(<r~ (#o)+) = N(cr~(6o)) = in that case. 
Either way, U is continuous on ]a~(9o), 9q[, hence constant by Theorem[TJ and 
so = L. 

Suppose that N(<r-(6 )) = c(<7 _ (0 o )) = c_. Then 

\N--c^\ = \N(9 Q -)-N(a-(9 Q ))\= Urn /' L(cj))d(j) (197) 



(0o) 



<(9 Q -<J-(9 ))\u\ m ^. (198) 

Recall 



\N_ 



^i + ^V v " V 2 7 



iSl (./nr. : 1^,1^1 



y/l + z- Wl + z-L— (199) 



7 + 1 Z n 



-max 



since that function is concave for z > and at z = 0. So, 

{9 Q -a-{9 Q ))>5'[j{U-9 Q ). (202) 

for some d'[ > 0. 

Suppose instead that N(a~(9o)) = 0. Then, 

AT_ = N(6q—) — N(a~(9o)) = lim / L{(j>)d<t> (203) 

»7/**o Ja-(e ) 

<{9 -cy-{9o))\u\ m ^. (204) 

From before, 



1 + z 



N- = c_ v V > CmiM /2-i. (205) 

Therefore, 

(So -o--(0 o )) > S > (206) 
in this case, for some S > 0. Then, taking 

S L =rmn(s , L ,d'l,-^-), (207) 

we obtain the desired result. (Recall that there is an upper limit to shock 
strength for the phase space under consideration). Analogous calculations yield 

the statement for backward-facing shocks. □ 
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9 Decomposition of the domain; regularity 

We now divide [0, 2ir) into the sets at which different behavior occurs. We shall 
assume that there are no stagnation points. 

C:=\e U is continuous at 6, \N(6)\ £ c{0) and ^ oj (208) 

Sf '■= 1 9 U has a forward facing shock at 6* j (209) 

Sb '■= |# U has a backward facing shock at o\ (210) 

S c := {0 | N{0) = 0} (211) 

Tl F := 1 6 U is continuous at 6, N(6) = c(0)} (212) 

K B := {0 U is continuous at 6, N{0) = -c{6)\ (213) 

Lemma 7. Recall Lemma\^ In a similar manner, TZp is uniformly separated 
from TZb U Sb U 5c, arad 72. b is uniformly separated from TZp U Sf U iSc- 

Proof. The proof is similar to the proof of Lemma [5] □ 

Lemma 8. There exists a constant Cs such that if 0q G Sf, and 
6<£\a-{6 ),a+(6 )[, then 

J(U; O ), \N(0 O +) - c(9 +)\, \N(8 -) - c(6 -)\ < C s \0 - Q \. (214) 
Similarly, if 0q G Sb, and £]cr~ (6q) , <j + {Qq)[, then 

J{U;6 O ),\N(0 O +) + c{0 o +)\,\N(0 o -) + c(0 Q -)\ < C s \0 - O \. (215) 
Proof. Suppose 0q <G Sf- From Theorem 01 we have 

S L J(U:0o) <mm{\a + (0o)-0ol\O o -a-(0o)\) < \0 - 6 \. (216) 
Therefore, 

J(U;6 )<5l 1 \6-e \. (217) 
From previous calculations, we have that 



1 + z 1 -^- - lj < c max + z^±± - 1 



The function on the right is at z = 0. concave down, and has derivative at 

~(7- 
4f 



z = equal to fssdl+ii. Therefore 



\N+ -c+\< Cmax i 7 + l) z< C" S J(U; Q ) < CsSl 1 ^ - 6 \. (219) 
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Also recall that 



\N- 



+ \i + ^V V 2 7 

< c max (vi + z - + . (221) 

This function is at z = 0, concave down, and has derivative at z = equal to 
Cmax ( I - ^r). Therefore, 

|JV- -c-l <ew^tl«<C^J(^;fl ) < Ps^l* - *o|. (222) 

Taking 



C 5 = H^.^A 1 .^*! 1 ) ( 223 ) 
gives the desired result. A similar argument works for # G <5s. □ 

Lemma 9. Recall that Sc is a finite set of points, and therefore can have no 
limit points. If 8 is a limit point of Sf, then 9 £ TZp- If 9 is a limit point of 
Sb, then 9 £ TLB- 
Proof. Consider {9 n } \ 9 be a strictly decreasing sequence in Sf- If 9 is a 
limit point of Sf, then 9 f\a~ (9 n ), cr + (9 n )[ for all n. For all n, choose 9' n £ 
}o-{6 n ),9 n [. Then, 

W(9' n ) - c(6'J\ = \N{9' n ) - c(0„-)|, since U is constant on ]<r~ (0 n ), n [, 

= \N(9' n ) - N(0 n -)\ + \N(0 n -) - c(9 n -)\ (224) 
< Mm^K - 8n\ + C s \9 - 9 n \, from Lemma E (225) 

= o{\e n -e\). (226) 

Thus we have a sequence converging to 9 such that 

Hm \N(6' n )-c(9' n )\=0. (227) 

n— too 

This eliminates the possibility of a shock or contact occurring at 9, and so U 
is continuous at 9. Therefore, 9 £ TZp by definition. A similar argument works 
for limit points of Sb- d 

Lemma 10. If 9 £ C, then U is constant on a neighborhood ]k~(9), k + (9)[ con- 
taining 9. (k {&) are taken to be maximal so that each is in either ScSftSb^f, 
or TZb-) 

Proof. From Lemma [SJ 9 is not a limit point of Sc,Sp, or Sb- Therefore, U is 
continuous on a neighborhood of 9. Therefore, if |iV(0)| 7^ c(9) or 0, then this 
will also be true on a neighborhood of 9. Then, Theorem [T] applies and so U 
is constant on some neighborhood containing 9. Clearly k ± (0) can be taken to 
satisfy the requirement in the statement of the lemma. □ 
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Lemma 11. Assume there are no stagnation points. Define, for 9 G [0, 2tt[, 

U S {6)= {U{4>+)-U{4>-)). (228) 

</>e[o,e[n(5 c u5FU5i3) 

Then Us is a right- continuous saltus function (so, by definition, is of bounded 
variation) . 

Proof. We have 

J2 \U{<t>+)-U{<t>-)\= J2 J (^)+ E \U(<t>+)-U(j>-)\ 

<j>e(S c US F US B ) cf>£(S F US B ) 065 c 

(229) 

<(2<5 i )- 1 E c < 

<j>eS F us B (fieSc 

(230) 

since the number of contacts is finite, the neighborhoods }a~ (cj)) , a + ((/>)[ are 
pairwisc disjoint, the phase space is compact, and the domain is compact. The 
BV norm only depends on the equation of state and the bounds for the phase 
space. □ 

Lemma 12. For any 9a € TZf U TZb, Us satisfies a Lipschitz estimate based at 
9q for 9 sufficiently close to 9q . That is, 

\Us(6)-Us{6o)\<Cs\6 -9 \, V '9 sufficiently close to 6 . (231) 

Moreover, the Lipschitz constant Cs is uniform in 9 and independent of U . 

Proof. Consider 8 € TZp, and 8 > 9 Q . We only need to consider other forward 
facing shocks occurring between 9 and 0, since TZp is uniformly separated from 
Sb and Sc- Suppose 9 ^]a~ ((f>) , o~ + (4>)[ for any G Sf- Then, 

\Us(0)-U s (9o)\< ]T J{U;d>) (232) 
e «t><e 

<(25 L )- 1 £ (a+(cb)-a-(cb))<(2S L )- 1 \8-9 l (233) 

0o«t><0 

since the }a~ (</>), a + (<p)[ are pairwise disjoint and contained in [9q,9] by as- 
sumption. If 9 £]a~ ((/>), (j>[ for some 4> S Sf, then the previous estimate holds 
for 9 = o~~(4>), and Us is constant on ]a~ (</>), <j>\ and so the result follows. If 
9 £ [<fi, o + (tj))[ for some <f> G Sf, then apply the previous estimate for 9 = a~{4>), 
and then 

\U S (8) ~ Us(8 )\ < (25 L )- 1 \(<T-(<t>)-9o) + J(U; ( f>) (234) 
< {{25L)- 1 +C s )\9-9 a \. (235) 

from Lemma [8] Take Cs = (25l)~ 1 + Cs- Similar arguments work for 9 < 9q, 
and 6» € TZ B - □ 
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Lemma 13. For every 9q £ TZf, there exists a neighborhood containing 9q such 
that 

u{6) sin 9 - v{9) cos 9 - c{9) (236) 

satisfies a Lipschitz condition based at 9q for all 9 in this neighborhood. The 
Lipschitz constant is uniform for all such 9q and is independent of U . That is, 

u(9) sin 6*o — v(9) cos6> — c(9) < M\9 — 9 \, V# sufficiently close to 9 , 

(237) 

(recalling that (N(6 ) - c(9 )) = for 9 Q g TZ F ). 
We have the similar estimate for 

u{9) sin 9 -v{8) cos 9 + c(9) (238) 

with 9q € TZb- 

Proof. Suppose #o € TZf- We first prove the desired estimate for \N(6) — c(9)\. 
Take 9 > 9q sufficiently close to 9q- This means that either 9 E C, with 
K~{9) e K F U S F , 9 e S F , or 9 e K F (from LcmmaEJ). 

Suppose 9 £ TZ F - 

Then \N(9)-c(9)\ = 0. 

Suppose 9 G Sf- 
Then, recalling 

cr" (6) = inf \n < N(4>) -c((j)) < c(</>), e]^ 6>[| , (239) 

it is clear that #o ^] f7_ (^)^[- Therefore, Lemma [5] applies and 

|JV(0) - c{9)\ = \N(9+) - c(9+)\ < C s \9 - 8 Q \, (240) 
(recalling we have made U right continuous everywhere). 

Finally, suppose 9 € C. If K~(9) G 7?-f, then 

|JV(e)-c(fl)| = |iV(e)- C («-(fl))| (241) 

< |JV(0) - 2V(«-(©))| + \N(k-(9)) - c( K -(9))\ (242) 
<Nnu*|K"(0)-0|+O (243) 
<Nm«|fl-flo|. (244) 

If «T(0) 6 5 F , then 

|JV(ff)-c(fl)| = |JV(d)-c(«-(fl))| (245) 

= |JV(fl)-JV(«-(fl))| + |JV(«-(ff))-c(«-(fl))| (246) 

< Nmax|0-K" (6)\+C s \9-6 \ (247) 
< N m ax|0-0o|+C s |0-0 o |- (248) 
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Taking M' = max(C<;, 



gives the desired estimate for \N(9) — c(9)\. Then, 



|u(0)sin0 o - v{9)cos6 - c(0)| = \N(9) - c(0)| + |u(0)(sin0 o -sin0) -ij(0)(cos0 o 

(249) 

<M'|0-0o| + 2|u| raax |0-0 o |. (250) 

Taking AI := M' + 2|u| max gives the desired result for > 0o, sufficiently close 
to O € 1l F . 



Similar arguments work for < 0q, and for 0q £ TZb- 



□ 



Wc now recall the concept of genuine nonlinearity. For fixed 0, the quantities 
N ± c are genuinely nonlinear in the sense that 



(N±c) u r ± (U;9) ^0 



(251) 



for all U. For simplicity we can compute the derivative in terms of the primitive 
variables, and recalling c = W7p we obtain 



(N±c) u r ± (U;l 



( Thf^ 



N±Ji-) Vu^iU-,9) 



sin 



P 



cos ± 



(252) 



1 Ax ^ 

2 V PP J 



2 p 



2 p 



COS [ 



2p J 



1 

U 
V 



p 





o 



i 



cos0 ±i£ ) 



V (7-1) 

/ 1 

±^sin( 

p 

=F- cos( 

" 2 

cr 



u +v 



1 c c 1 C J 

T 2p p 2~p 

1 c c 7 c 
T 2p p 2p 

±^(7 + l) £ t^O. 

2 p 



u ± c sin 
u =F ccos 
\ H±Nc J 







1 



-(7 - l)u -(7 - l)n 7 - 1 / 



uic sin ^ 
u =p c COS 
\ H±Nc J 



(253) 
(254) 
(255) 



The Roe linearization for the full polytropic Euler equations has the advan- 
tage that it is simply the matrices f§ and ffj evaluated at some appropriately 
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averaged state U. It takes the form (see |T09j ) 
A(U.,U + ;9) :=A(U;9) = 
( 



N ■ 



sin 9^- (u 2 + v 2 ) -TiN_ 

■ cos 9 :L ^-(u 2 + v 2 ) — UN TJsin0 + (7 - l)ucos( 



sin 9 
(2 — 7)usin( 



V (^{u 2 +v 2 )-h) N 
We define 



hsm9 - (7- l)Nu 



(256) 

— COS0 

-wcos0 — (7 — l)zJsin0 

N + (7- 2)?Jcos0 
-hcosO- (7 - 1)A% 



° \ 

sin #(7 - 1) 

- cos 0(7 — 1) 
7 iV / 



1 



- = — L ^~ + ^ 2 + v 2 ). 
P 7- 1 P 2 



to be the total specific enthalpy per unit mass, so that 

c 2 = (7 - 1) I 



is the sound speed. The averaged quantities are defined as 



VP- + VP+ 

V = / 1 r— ' 

N = u sin — v cos 9, 



—2 1 —2 

u + v 



VP- + VP+ 

(7-1) (h- 



The eigenvalues of A(U; 9) are 



N±c 1 N, N, 



(257) 
(258) 

(259) 

(260) 
(261) 
(262) 

(263) 
(264) 



and it has a full basis of eigenvectors. Moreover, it is clear that A(U; 9) is 
a smooth function of C/±, the eigenvalues are smooth functions of U±, (away 
from p = of course), and by direct inspection of the eigenvectors (they are 
not needed here, but expressions for them are available) they too are smooth 
functions of U± . Define the left and right eigenvectors so that 



A(U; 0)^(17 
^(77; 9)A(U 
A(U; 9)r l (U 
9)A(U 
l a {U; 9)r a {U 



9) = (N±c)r ± (U; 0), 

0) = (N±c)l ± {U;9), 

9) = Nr l (U;9) for £ = 1,2, 

0) = Nl l (U;9), for £ = 1,2. 

9) = S aP , a,0 = +,-,l,2. 



(265) 
(266) 
(267) 
(268) 
(269) 
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Theorem 5. For any 9q G IZfUIZb, there is a neighborhood containing 9q such 
that U satisfies a Lipschitz condition based at 9q for all 9 in this neighborhood. 
The Lipschitz constant is uniform for all such 9q and is independent ofU. That 
is, 



U{9) - U(0 O ) < M'\9 - O |,V0 sufficiently close to 9 - 



Proof. Suppose 9q £ Tip. 
Recall that 



A(U(0 o ), U(9); 9 ) (U(0) - U(9 )) = 0(\9 - 9 \). 



(270) 



(271) 



Denote l + (U(9o),U(9);9o) = l + (U;9 ) where the average is taken between 
U(9 ) and £7(0). Left multiply by l + (U{9 ), U{9); O ) to obtain 



{N + c)l+(U(6 ),U(6);6 Q ) (U(9) - U(9 )) = 0(\9 - 9 Q \), 



(272) 



where the averages are taken between U(8o) and U{9). Since 6> £ TZf, U is 
continuous at 9q. Therefore for 9 sufficiently close to 9$, N + c is uniformly 
bounded away from zero, since N — c will be approaching zero. Therefore, 



i+(u(e ),u(8 1 y,e )(u(9)-u(e )), 



(273) 



satisfies a Lipschitz estimate based at 9q for 9 sufficiently close to 9q, with 
Lipschitz constant uniformly bounded above by the bounds on the phase space, 
and proportional c~ s . Similarly, 



9 i ^ l l (U(9 ), U(9); 9 ) (U(9) - U(9 Q )) , i = 1, 2, 



(274) 



also satisfies a Lipschitz estimate based at 9 for 9 sufficiently close to 9 , with 
a similar upper bound on the Lipschitz constant. We claim that 



W^g{W) 



9 2 {W) 


( 




g 3 (W) 




V g\w) ) 


\ 



w 



sin 9q 



cos 9 



z{W) \ 



1+{U(8 ) 7 W;9 )(W - U{8 )) 
l 1 (U(8 ),W;8 )(W ~U(9 ) 



W-U(9 )) ) 



(275) 



defines a diffeomorphism for W sufficiently close to U(8q), when combined with 
the previous lemma will prove the claim. Notice that 



9w(U(0 o )) = 



/ (iisin6> - vcos9 - c)u\ v , e s \ 
l+(U(9 Q );9 ) 
l 1 (U(9o);9 ) 
V 1 2 (U(9 );9 ) ) 



(276) 



Then, if 



g l w (U(9 ))z = 0,i = 2,3,4, 



(277) 
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this implies z \\ r (U(9o);9 ). But since 

((usin0 o -ucos6> -c^r"^^)) ^ 0, (278) 

V / U(8 ) 

(by genuine nonlinearity) this implies z = 0. Therefore, for W sufficiently 
close to U(9o), g is a diffeomorphism. Since U(9) approaches U(8 ) as 9 — >• O 
by continuity, and g(U(6)) satisfies the Lipschitz estimate based at 9q for 9 
sufficiently close to 9q, U itself must satisfy a Lipschitz estimate based at 9q 
for 9 sufficiently close to 9$. The C 1 norm of g is bounded uniformly above by 
the phase space bounds and equation of state, as is the Lipschitz constant for 
g(U(6)), and so the Lipschitz constant for U is as well. □ 

Theorem 6. Assuming that there are no stagnation points, and that density 
and internal energy remain bounded away from zero, any L°° weak, steady, self 
similar solution to the 2-d full polytropic Euler equations must be of bounded 
variation. Moreover, U can be decomposed as 

U = U L + U Sl (279) 

where Ul is Lipschitz with constant independent of U , and Us is a saltus func- 
tion of bounded variation, with total variation independent of U . (Note these 
constants will depend on the equation of state, the lower bound on density and 
internal energy, and the L°° norm of U .) Note that this implies U is a special 
function of bounded variation, since the Cantor part vanishes. Moreover, the 
absolutely continuous part is in fact Lipschitz. 

Proof. The statement about Us has been covered in previous lemmas. We claim 
that for any 9o, Ul satisfies a Lipschitz estimate based at 9o for 9 sufficiently 
close to £>o- 



If #o G C, then U is constant on a neighborhood containing 9q, and since there 
are no shocks or contacts it is clear that Ul '■— U — Us is constant and thus 
satisfies a Lipschitz estimate based at 9$ with constant 0. 

If $o G Sf U Sb U iSc, then the jump at #o is accounted for in Us, and so 
Ul is constant on some neighborhood containing #0j and so satisfies a Lipschitz 
estimate based at 9q with constant 0. 

If 6*o G TZf U TZb, then, for 9 sufficiently close to 6q, we have from Lemma 
[Hand Theorem [5] that 

\U L {9) - U L (0 O )\ < \U(9) - U(6 )\ + \U S (9) - U s (9 )\ (280) 
< C s \9 - O I + M'\6 - O | := C L \0 - 9 \. (281) 

Pick any 9\ and 9i. Consider the open cover 

U (282) 

4>e[e u 6a] 



33 



where for any 9, H(9) is the neighborhood for which we have a Lipschitz estimate 
based at 9 with Lipschitz constant uniformly bounded by Cl- This has a finite 
subcover 



/v 



(J (283) 

n=l 

Then, adding in Q(6i) and Q(9 2 ), we can express 

\U(9 2 ) - U(9 1 )\ < \U( m ) - U(9 1 )\ + \U(fa) - U(m)\ + Mm) - U(fa)\ (284) 
\U(fa) - U{ m )\ + ... + \U(9 2 ) - U( m+1 )\ (285) 
<C L \9 2 -9 1 \, (286) 

where fa-i < r\i < fa, and ry, € Sl(fa-i) n Q(fa) for i = 2, .., AT. (We take 6*i < 
?7i < 01,17], G Q(6»i)nQ(0i), and fay < i] N+1 < 9 2 , with 77^+1 G n{<j) N )r\fl(9 2 ).) 
Thus Ul is Lipschitz on all of [0, 2iv), and the rest follows. □ 



10 Structure of flows 

We now prove some results about the structure of possible solutions, and present 
several examples. 

Wc begin by decomposing the domain into a finite number of sectors. Denote 
the points in Sc as 9\, 9 2 , 9jq so that 



< 6»i < 9 2 < ... < 9 N < 2tt. (287) 

(Recall that 9 £ Sc means that N(9) = 0, and Sc is a finite set by Lemma 
2J) Define the sectors lj, for i = 1, .., N, as 

li:=[0i,0»+i]> (288) 
(taking 9^ + \ = 9\ to unify the notation). 

We say that Ii is a, forward sector if N\i i > 0, and that is a backward sector 
if ^l/i < 0. By construction, each sector will either be forward or backward, 
and N will be positive on the interior of a forward sector, and negative on the 
interior of a backward sector. Moreover, L is continuous on the interior of each 
sector, since L is continuous at shocks and there are no contacts in the interior 
of a sector. 

For all the figures in the remainder of the paper, the flow direction is what 
is indicated. The length of the arrows is not meant to suggest anything about 
the length of the velocity vectors. 
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Figure 1: In a forward sector 7j = L(9i+) > and L{Q i+1 — ) < 0. L 

is monotonically decreasing on ]#.;, and equal to zero at a unique (9;. The 

flow is constant on ]9i, cr + (6i)[ and ]a~ [. If 7^ 0j + 7r, then there 

must be some wave structure in the grey shaded region. 

Lemma 14. (See Figure\J\) Suppose Ii is a forward sector. Then, L(6i+) > 0. 
L(9i+i—) < 0, and L is monotonically decreasing on ]9i,6i+±[. Similarly, if Ii 
is a backward sector then L(9i+) < 0, L[9i-^i— ) > 0, and L is monotonically 
increasing on ]6i,6i+i[. Therefore, there exists a unigue 8i G.]6i,9i+i[ such that 
L(?i) = 0. 

Proof. Consider the strong form of the conservation of mass and tangential 
momentum equations, 

(pN) e = pL, (289) 
(pLN)e = pi? - pN 2 . (290) 

Manipulating these, we obtain for any point of differentiability on the interior 
of Ii that 

LgpN + L(pN) e =pL 2 -pN 2 , (291) 
L e pN + pi? = pL 2 - pN 2 (292) 
U = -N, (293) 

since p is bounded away from zero by assumption, and N ^ on the interior of 
Ii. 

Theorem |6] shows that U is Lipschitz almost everywhere (since Us is con- 
stant except on at most a countable, discrete set), and since jumps in Us on 
the interior of Ii must be shocks (not contacts) , Ls is constant on the interior 
of Ii. Therefore, L is Lipschitz (hence differentiable almost everywhere) on the 
interior of It, and so the fundamental theorem of calculus can be applied to 
L. Therefore (|293|) shows that L is monotonically decreasing (increasing) on 
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if Ii is a forward (backward) sector. 

Recall that when U is constant, Ng = L. Also recall that by Theorem [3] there 
exist a + (0i) > Oi and a~(6i + i) < 0i + \ such that U is constant on ]0i,a + (9i)[ 
and on ]a~ Therefore, the following right limits arc defined and we 
have for small 5 > that 

N(0 i + S)-N[0 i +)=N(0 i + 6) = / (294) 

Since U is constant on ]9i,6i + <5[ and there are no stagnation points ^ 
\&(0i+)\ = \N(6i+)\ + \L(6i+)\ = \L(6i+)\. Therefore, by continuity, sgn(L) 
is constant on ]6i,6i + S[, and so sgn (N(9i + <5)) = sgn (L(6i+)\. Therefore, 
for a forward sector, L(6i+) > 0. Similar arguments work for L{6i+i— ) and for 
backward sectors. 

Since L(6i+) and L(6%+i—) must have opposite signs, and L is monotonic on 
there is a unique Oi such that L(6i) = 0. □ 

Now we define a Prandtl- Meyer wave. A forward Prandtl- Meyer wave is 
a closed interval [a,/?] such that N(8) = c(9) for all £ [a,/3]. A backward 
Prandtl-Meyer wave is the same except that N{&) = —c(ff). Moreover, U is 
diffcrcntiable almost everywhere on ]«,/?[, and Ug is in the kernel of A(U;0) 
when it is defined. This follows from the the strong form of the Euler equations 

A(U;0)U =O. (295) 

It is well known that p = A{S)p 1 for a polytropic gas, where A(S) is a func- 
tion of the entropy, and that away from discontinuities S is constant. Therefore, 
for 8 in the interior of a forward Prandtl-meyer wave we have that 

(pN) e ^(pc)g = ^A(Sh{p^ 1 ) g (296) 

= V^h^P^Pe - PL, (297) 

and so 

sgn(p e ) = Bgn(L). (298) 

Therefore, for a forward Prandtl-Meyer wave, as the gas particles pass through 
it (corresponding to the decreasing direction by our choice of coordinates), 
p increases if L is negative, and decreases if L is positive. Since the flow is 
isentropic, 

sgn(pe) = sgn(p e ) = sgn(L). (299) 

Similar calculations for backward waves can be done. Therefore, in light of 
Lemma H~4l we have the classifications: 
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• forward expansion wave: a < (3 < Oi, N(0) = c(8),L(0) > for all 
0e [a, (3], 

• forward compression wave: 9i < a < f3, N(8) = c(9),L(8) < for all 
9€[a,/3], 

• backward expansion wave: 9i < a < /?, N(9) = —c(9),L(9) > for all 
e [a,0\, 

• backward compression wave: a < (3 < 9i, N(8) = —c(8),L(9) < for all 
9€ [a, 13}. 

It is possible to join a forward compression wave between [ol\, fti] to a forward 
expansion wave between [ao,/3o] if Pa = 9i = ol\. In that case the compres- 
sion wave ends when the flow is precisely sonic at 8i, and an expansion wave 
immediately starts at 8i. We have the following theorem. 















L < 0\ 






wy l > o 






— ► \^ 






9i+i 











Figure 2: In a forward sector Zj = [8i,8i + i], L > on and L < on 

]8i,8 i+ i[. The flow is constant on }9i,cr + (8i)[ and }<r~ (Oi+i), @i+i [• There is at 
most one shock or rarefaction in [<r + (9i), and possibly infinitely many shocks 
and compression waves in a~ {9i + i)\. However, there cannot be consecutive 
compression waves. In this particular example, the flow consists of a compression 
wave and two shocks in the L < part, and a rarefaction wave in the L > 
part. 

Theorem 7. ( See Figure [H ) Suppose Ii is a forward sector, and that U is 
continuous on an open interval B c]0j,#i+i[. (In this case L < on B.) Then, 
either U is constant on this open interval or constant on either side of a single 
forward compression wave. 

Similarly, in a backward sector, on an open interval B £]6i, 9i[ on which U is 
continuous, U must be constant or constant on either side of a single backward 
compression wave. 



37 



Proof. Suppose Ii is a forward sector. Since U is continuous on B, the set on 
which N(9) = c{9) is closed in B. Therefore, its complement in B is a countable 
union of open intervals, on which U is constant by Theorem Q] (since N{ff) can- 
not be or —c(9) in the interior of Ii). Since L(9) is negative on B, and Ng = L 
on this complement, N(6) = c(9) can be satisfied at at most one endpoint of an 
open interval in the complement. Therefore at least one endpoint must be an 
endpoint of B, making TZf D B a closed interval in B. 

Therocm [5] shows that U is Lipschitz on B, since B C C U TZf and thus Us 
is constant on B. Therefore, U is differentiable almost everywhere in TZf H £? 
and the strong form of the equations implies Ug is in the kernel of A(U; 9) 
everywhere it is defined, and so TZf H -B defines a forward Prandtl-Meyer wave. 
The fact that B C]Oi,0i+i[ shows it must be a forward compression wave. The 
argument for a backward sector is similar. □ 

Note that there may be multiple forward compression waves in a forward 
sector - this theorem requires only that there is at least one forward facing 
shock in between. Since L is negative on on any interval on which U 

is constant TV is decreasing. In a forward sector, this corresponds to TV increasing 
along particle paths of the gas particles. Therefore, upon exiting a compression 
wave, the normal velocity is sonic, but as the gas particles continue traveling in 
the negative 6 direction, the normal velocity increases and becomes supersonic, 
leading to the possibility of a forward facing shock, which upon exit the normal 
velocity will be subsonic. Normal velocity can then increase along particle paths 
back to the sound speed, and the gas can enter another compression wave. 

Theorem 8. (See Figured) Suppose Ii is a forward sector. Then, on ]6i,0i], 
exactly one of the following is true: 

• U is constant on either side of a forward facing shock, 

• U is constant on either side of a forward expansion wave, 

• U has an expansion wave on [a, and is constant on ]9i,a[, 

• U has a normal shock (that is, I = 0) at 9i and is constant on]9i,8i[ , 

• U is constant on ]9i,9i]. 

We have the similar statement if Ii is a backward sector, for the interval [9i, 9i + \ [. 

Proof. Suppose U has a shock at 9q G]9i,9i\. Then, we know that N(9q — ) < 
c(9q—). Recall that there exists a~(9o) < 9q such that U is constant on 
}o--(6 ),6 [, and that N(a~(9 ) + ) = or c(<7 - (0 o ) + ). However, Ng = L 
on ](7~ (6*o), 9[, and on this interval I > 0, and so TV decreases as 9 decreases, 
and so TV^^^o) + ) = 0, making a~(9o) = 9\. Therefore, there can be no 
shocks in ]9i,9o[, and U is constant on ]0»,0o[ by Theorem [TJ Similar argu- 
ments show that a + (9o) > 9i (since for TV to be sonic it must decrease from 
N(9q) > c(9q), which is impossible since Ng = I > on ]9q, a + (9q)[) and so 
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either the first or fourth statement is true. 

If there is not a shock, then U is continuous on this interval, and similar 
arguments as in the proof of Theorem [7] show that there can be at most one 
expansion wave, and we are done. Similar arguments work in backward sectors. 

□ 

Examples with infinitely many shocks can be constructed (these theorems 
show that they must occur in the parts of the sectors where L < 0), or with 
infinitely many shocks interspersed with compression waves (with the restriction 
that compression waves cannot occur consecutively, by Theorem [7]). Therefore, 
since infinitely many discontinuities may occur BV is the sharpest commonly 
used function space we may use. 



11 Maximum number of contacts 

We note for a both rarefaction and compression waves that the velocity turns 
towards the origin as the gas particles travel through the wave. This can be 
seen by manipulating 

( P N)g = pL, (300) 
( P N 2 +p) e = 2 P NL, (301) 

to obtain the following. 

(pN)gN + (pN)Ng +pg = 2pNL, (302) 
pLN + (pN)Ng + Pe = 2pNL, (303) 

NN e - LN = -—, (304) 
P 

LN — NN g = — . (305) 
P 

We consider the angle of the flow, 4> := Z(u, v), as in Figure [3] , as in |LL87j . 

Since we are considering N > 0, Z(N,L) € (— tt/2, 7r/2), and so Z(N,L) = 
arctan (-4 ) . Then, 

Z(u, v) = 6 - | + arctan (J^J (306) 
= - arctan ( — J . (307) 
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Then, 

, , 1 LNg-NLg 

^^TTTZ^ — v — 



(308) 



LNg + N 2 
L 2 + N 
L 2 - LNg 
~ L 2 + N 2 
_ L(L - Ng) 
L 2 +N 2 

_ Lpg 

~ pN(N 2 + L 2 ) ' 
However, recall that in a forward wave 



' (309) 



(310) 
(311) 
(312) 



sgn(pe) = sgn(L), (313) 

and so 4>g is positive. Therefore, as the gas particles travel through the shock, 
9 decreases, and so Z.(u,v) decreases as well. 

For backward waves, 



<f>=9- arctan ( — J + tt, (314) 

giving the same expression for cj>g. But in this case 

sgn(j%) = -sgn(L), (315) 

but N = —c and so <pg is still positive. However, for backward waves the gas 
particles move in the increasing direction, and so the flow still turns towards 
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the origin as the gas particles travel through the wave. See Figure 2] for some 
examples of waves. 

For shocks, since L is continuous and \N\ decreases as the gas particles pass 
through the shocks, the flow is turned away from the origin if L is positive, and 
toward the origin if L is negative. 



^\~~^\ \ 


L < 

•• L > 


0i " 


>- — 1 




L < 

e 


L > O^v 

1 



Figure 4: An example with two sectors. I\ = [61,82] is on the bottom and is a 
backward sector, I2 — [62, #1] is on the top and is a forward sector. In I\, there 
is a single shock in the L > part. In I2, the gas passes through a shock, a 
compression wave, then a shock in the L < part, and a shock in the L > 
part. As the gas particles travel through shocks, the flow is turned toward the 
shock line, and compression waves turn the flow toward the origin. Rarefaction 
waves would turn the flow toward the origin as well. 

Theorem 9. For 7 = 1.4, there can be a maximum of two sectors. For other 
values of '7 > 1, there can be up to three sectors, but there are no values 0/7 > 1 
that lead to flows with more than three sectors. 

Proof. Choose coordinates so that there is a forward sector / = [a, it] where 
< a < 7T, so that N(n) = N(a) = 0. Then, by Lemma [H L(a+) > 0, 
L(ir—) < 0, and so the flow needs to be turned away from the origin by an 
angle of a. Following Qj] denote 8 the unique value between a and it such that 
L{8) = 0, and recall that L is positive on ]a, 8[, and negative on ]8, w[. 

By Theorem [71 the discussion after it, and the discussion preceding this 
theorem, any shocks or compression waves on ]8, n[ turn the flow towards the 
origin. Therefore, Z(u(8),v(8)^j < 0. Since we are interested in finding the 
maximum possible a that the flow can be turned upwards, the best possible sit- 
uation is for there to be no compression waves or shocks on ]8, tt[, which yields 
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6= f,Z(u(7r/2),t;(7r/2)) =0. 

If there is a shock at ^ , then it is a normal shock and so the flow is absolutely 
subsonic, and thus constant, for ]a, ^ [. Therefore a — since the flow can never 
be turned away from the origin. 

Therefore, to accomplish the maximum upwards turning, the flow should be 
constant on [§, 7r[- By Theorem[51 the flow is either constant (again resulting in 
a = 0), has exactly one rarefaction, or exactly one shock on ]a, -|[. A rarefaction 
wave turns the flow towards the origin, resulting in a < 0, and so there must 
be a single shock to accomplish a > 0. 

Using the well known 8 — (3 — M equation (see |K05| , Chapter 4) to relate 
the incident Mach number M := Mt ; the turning angle a, and the shock angle 
8, we have that 

/ 2 cot 6(M 2 sin 2 9 - 1)\ , . 

a = arctan ; r „ , — . (316) 

V Af 2 ( 7 + cos(26<)) + 2 J y ' 

It is well known (see |CF48j ) that the curves a(6) for fixed values of M all lie 
below the limiting case M oo, and solving for the maximum a yields 

aw* = arcsin (j^j ■ ( 317 ) 

The flow can only be turned upward when L > 0, so there can never be more 
than three sectors since a max = f is only attained in the limit as 7 \ 1. For 
7 = 1.4, a max w 45.5°, and so the flow cannot turn the required 60° needed to 
have more than two sectors. For 1 < 7 < 1.15, a max > 60°, and so there will 
exist finite incoming Mach numbers for which the flow can turn 60°, allowing 
for the existence flows with three sectors for some values of 7 > 1 . □ 
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Figure 5: For 7 < 1.15, the maximum turning angle is greater than 60°. There- 
fore there exist flows with three contact discontinuities, such as the one above. 
In this example each sector has one shock in the region where L > 0, causing 
the flow to turn away from the origin. 
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